Self-cooling of a movable mirror to the ground state using radiation pressure 
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We show that one can cool a micro-mechanical oscillator to its quantum ground state using 
radiation pressure in an appropriately detuned cavity (self-cooling). From a simple theory based on 
Heisenberg-Langevin equations we find that optimal self-cooling occurs in the good cavity regime, 
when the cavity bandwidth is smaller than the mechanical frequency, but still larger than the 
effective mechanical damping. In this case the intracavity field and the vibrational mechanical 
mode coherently exchange their fluctuations. We also present dynamical calculations which show 
how to access the mirror final temperature from the fluctuations of the field reflected by the cavity. 
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Cooling of mechanical resonators both at the micro- 
and at the macro-level, is an important technical chal- 
lenge in various fields of physics, such as ultra-high pre- 
cision measurements [H, and detection of gravitational 
waves 0. It is also a prerequisite for any possible use 
of optomechanical systems for quantum information pro- 
cessmg 0,13. Active noise control techniques have been 
proposed to reduce their thermal noise and bring the os- 
cillator motion to its ground state Recently, var- 
ious experiments have demonstrated significant cooling 
of the vibrational mode coupled to an optical cavity 
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The experiments of 
12l . |13| in particular have adopted the so- 



Refs. 

called back-action 1^, or self-cooling scheme, in which 
the radiation pressure of an appropriately detuned cav- 
ity interacting with the mechanical oscillator, is used. In 
all these experiments, however, the resulting equilibrium 
state of the oscillator is classical, because the new mean 
excitation number is still much larger than one. There- 
fore it is important to establish the fundamental limits 
of self-cooling and if it is able to cool a mechanical de- 
gree of freedom down to its quantum ground state. We 
will describe the system dynamics in terms of quantum 
Langevin equations (QLE) for the Heisenberg operators 
of the system and we will show that ground state self- 
cooling is possible in the good cavity regime, i.e. when 
the mechanical resonance frequency is larger than the 
cavity bandwidth, and provided the cavity detuning and 
bandwidth are appropriately chosen so as to maximize 
the scattering of noisy photons into the cavity mode. 

We consider a driven optical cavity coupled by radia- 
tion pressure to a micromechanical oscillator. The typ- 
ical experimental configuration is a Fabry-Perot cavity 
with one mirror much lighter than the other (see e.g. 
0) B B HO) [iH EBl), but our treatment apphes to other 
optomechanical systems, such as the silica toroidal micro- 
cavity of Refs. [13, 113 ■ Radiation pressure typically ex- 
cites several mechanical degrees of freedom of the system 



with different resonant frequencies. However, we consider 
the detection in a narrow frequency bandwidth, contain- 
ing a single mechanical resonance peak. In this case the 
oscillator dynamics can be well approximated by that of 
a single harmonic oscillator with dimensionless position 
q{t), and momentum p{t), {[q,p] = 2i), frequency rim, 
mass M, and damping F. The coupled dynamics of the 
mirror with the intracavity field mode a (t) in the frame 
rotating at the frequency of the driving laser uj^, is de- 
scribed by the QLE 



p = —i}„iq^Tp + GV2a^a + ^ 
a — ~ {k -\- i/!^c) o, + iGaq/V2 + 
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where G = {ujc/ L) ^ h/AdVlm is the opto-mechanical 
coupling constant, uJc and L being the cavity resonance 
frequency and length, respectively, a™ (t) is the input 
field, satisfying (a'"(i)a'"t(t')) - |a'"|^ = S{t~t'), with 
a"^ the incident mean field, k is the damping rate of 
the cavity mode, A^, — lUc ~ ujl the cavity detuning and 
^ the noise operator accounting for the mirror Brown- 
ian motion at thermal equilibrium at temperature T . 

Steady state analysis - In steady state the mean in- 
tracavity field a is given by a = •\/2Ka'"/(K -I- iA), 
where A is the mean detuning of the cavity, given by 
A = Ac — Ani = Ac — G^|ap/r2„i. These two coupled 
equations give a third-order relation between a and A 
which leads to the well-known bistable behavior of a cav- 
ity with a movable mirror [l9| . The stability condition 
of the system can be written as -\- A'^ -\- 2AA„i > 0. 

If we then consider the linearized fiuctuations of the 
various operators around the steady state, we get lin- 
earized QLE which can be solved by moving to the fre- 
quency domain. The Fourier transform of the position 
fluctuations can then be simply expressed as the sum of 
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a thermal noise term and a radiation pressure term 



(4) 



the response to the noise terms being given by an effective 
susceptibihty 



(5) 



with xl^] = Mfil^ [1 - n'^/nl^ - irn/nf^] and D[n] = 

(1 — ift/n)'^ + (A/k)^. The thermal noise force is sinrgly 
Jt[^^] = Mr2mC[f^] and the radiation pressure force [20| 



with x"^ = a'" + a'"t and y'" = i(a'"t - a'"). An exact 
expression for the mirror variances is then obtained by 
integrating the contributions of these two forces to the 
noise spectrum (cj = fl/rim) 



(6) 



where we have used the fact thatp 
the noise spectrum is given by 



q/^lm and where 
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with 6 = 1^ = A/k, ifni = A„i /k and Q = Slm/F. 

These two exact expressions for the variances coincide 
with Eq. 5 of Ref. [2l[. Cooling to the ground state 
implies reaching Aq^ — 1 simultaneously with Ap2 — i. 
Fig.[T]shows both variances as a function of b, for a cavity 
detuning equal to one bandwidth [ip = b), which shows 
that, starting from a mean thermal excitation number 
= {exp{hil„i/kBT} — equal to 102, it is possi- 
ble to achieve a lower than unity final excitation number 
state in the good cavity limit: 6 > 1, i.e. > 




FIG. 1: Normalized variances versus 6, for a detuning </? 
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' 0.15 for b . 
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0.1. The final excitation number 



Physical interpretation in the adiabatic limit - To get 
more insight into the cooling mechanism we will now 
derive an analytical expression for the frequency inte- 
grals of The latter may be considerably simpli- 



fied by noting that, for high-Q cavities and under cer- 
tain conditions that we will detail, the mirror response 
is peaked around the natural mechanical resonance fre- 
quency. One may then express the effective susceptibility 
as xP"^ - Mn'^[l-n'^/n'^-i/Q] with an effective qual- 
ity factor Q = and effective resonance frequency and 
relaxation rate 



n = rim [i - 2(/3<^„i Re[r>[ri™] 

f = T [1 + 2(f^nlQ MD[n,nr^ 
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(8) 
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The effect of the frequency-shift in the resonance is typ- 
ically negligible, so that one may assume f2 ~ Qm- The 
relaxation rate of the mirror is on the contrary strongly 
affected by the presence of the cavity field and, depend- 
ing on the sign of the cavity detuning, may either be 
enhanced or reduced, result ing in either cooling or heat- 
ing. As pointed out in [l^, [Hi , the most interesting 
regime for self-cooling is when the effective damping rate 
r is strongly enhanced, but still less than the cavity band- 
width and when the effective mechanical Q is still larger 
than one, F <C F ^ k. In this case, the effective sus- 
ceptibility is still peaked around uj = 1 and we can well 
approximate the smoothly varying function of lu inside 
the square brackets by taking its value at ~ 1. One 
then gets an analytical expression for the variance 
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Eq. (|T0|) is the basic equation for our analysis of the quan- 
tum limits of self-cooling, and it is possible to see that 
it coincides with the results of [13, [23|, w hich were ob- 
tained with a different method. Ref. I2ll| instead gives 
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slightly more general expressions for the two variances, 
which, however, reduce to that of Eq. pUj) when f2„ is 
large enough, i.e, flm > T. 

Depending on the respective frequencies involved, one 
can adopt different strategies to cool the mirror down to 
the ground state. In particular, the cavity detuning has 
to be well-chosen in order to optimize the self-cooling. 
Fig. [5] shows the variation of the variance with the cavity 
detuning, as calculated from the exact solution of Eq. ([7]) 
and from the harmonic approximation (|10p . The analyt- 
ical result is in good agreement with the numerical result 
and we will use it as a basis for the cooling optimization 
discussion. 




FIG. 2: Normalized variance versus cavity detuning [dashed: 
exact result of Eq. ([7]), plain: approximate result of Eq. (|10p ]. 
Parameters: Q = 10'', = 10^ 0.1, 6 = 10. 



In order to interpret Eq. fTU]) we rewrite it as 



Ag2 



1 + 2n\ 
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2ipb 



1 + ffnlQ 1 + f-i^nlQ 

or, equivalcntly, 

Ag2 = (1 - ,y) a4 + Aql. 



(11) 
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Aq'^ = 1-1- and Aq|, = (1 + + ip^)/2ipb represent 
the thermal noise-induced and the cavity field-induced 
contributions, respectively, to the position fluctuations 
of the mirror. 
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1 + f^nlQ 

is the weight between those two quantities, with 
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The physical interpretation of Eq. (fT^ is that the cav- 
ity field and the mirror vibrational mode coherently ex- 
change their fluctuations during the interaction, all the 
more so that the quantity rj is close to unity. The thermal 
excess noise is then replaced by the cavity field fluctua- 
tions, which are those of the vacuum, hence resulting in 



decreasing the effective temperature of the mirror. The 
quantity rj naturally appears as a quantum state trans- 
fer efficiency between the field and the vibrational mode, 
whereas the product ^PniQ can be interpreted as a co- 
herent coupling strength, quite similarly to the coopera- 
tivity parameter, which appears in quantum state trans- 
fer schemes between optical fields and atomic ensembles 
25, 26, 27]. Since the mirror is initially in a noisy ther- 
mal state, the cavity field acts as a thermal noise "eater" 
during the interaction, and allows to reach a much lower 
effective temperature for the mirror in steady state. This 
situation is also quite similar to cavity assisted Doppler- 
cooling of atoms [l!], for which the temperature of the 
atoms is decreased by enhancing the scattering of pho- 
tons into the cavity mode. 

If one starts with a very large thermal excess noise, 
Eq. (|lip clearly shows that the thermal noise contribution 
will be suppressed by a factor f(pniQ, which is typically 
large. Minimizing the mirror final temperature is then 
equivalent to maximizing the transfer efficiency rj, and, 
for a given value of fniQ, maximizing / with respect to 
if. This gives an optimal cavity detuning equal to 



= [(52 - 1 + 2Vl + &2 + 64)/3]^ 



(15) 



The most favorable conditions to reach the quantum 
ground state of the mechanical oscillator then occur for 
a cavity bandwidth smaller than the mechanical reso- 
nance. Indeed, assuming b ^ 1, one has ~ &, /* ~ 1 
and for a sufficiently high Q-factor {fniQ ^ 2n^ + 1), 
the resulting mirror fluctuations are then mostly given 
by the radiation pressure: 



(Aq^ 



(16) 



and the normalized variance tends to unity for & 3> 1. 
It is therefore possible to reach the mechanical ground 
state using self-cooling in the good cavity limit. 

However, there exists an optimal bandwidth that 
minimizes the mirror final temperature. Indeed, as 
one decreases the cavity bandwidth (for a fixed value 
of the mechanical resonance frequency), the effective 
response time of the mirror becomes comparable to 
that of the intracavity field and the previous adia- 
batic approximation is no longer valid. One can also 
show that this situation is detrimental to the noise 
transfer between the field and the mirror [2^, [2^. As 
mentioned previously and as can be seen in Fig. [1] the 
adiabatic approximation breaks down when ipni^m ^ 2k. 

Dynamical situation - We now address the dynamical 
behavior of the coupled field-mirror system. Under good 
self-cooling conditions one expect the mirror initial ther- 
mal noise to be mainly transferred to the cavity field 
during the cooling phase, and hence to the field leak- 
ing out of the cavity. Numerically solving the linearized 
Langevin equations in time allows to compute the two- 
time correlation functions oiq{t) and p{t) , as well as those 
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FIG. 3: (a) Temporal evolution of the mirror position and 
momentum normalized variances when the cooling field is 
on from t = {T~^ units), (b) Temporal evolution of the 
outgoing field variances Cx,y{t,t) (in k units). Parameters: 
K = lO^r, Q = 10*, nip = l6^ ifini = 0.1, V3 = fo = 10. 

of the outgoing field quadratures, a;°"* and y™*, where 
a°"* = \f2Ka — a™. To simplify, we assume a step-like in- 
jection of a coherent state cooling field into the cavity at 
i = and start with an occupation number = 10^ for 
the mirror [2^ . The resulting mirror variances at time 



t are represented in Fig. [3l as well as the temporal evo- 
lution of the equal-time variances of the outgoing field 
Cx.y(t, t), obtained from exact numerical calculations. It 
appears clearly that, during the cooling phase, thermal 
noise is taken from the mirror and transferred to the cav- 
ity field. The amount of thermal noise taken from the 
mirror can be deduced by measuring the outgoing field 
fluctuations. Indeed, for = A and in the adiabatic 
limit r <C r ^ K, it can be shown that the outgoing field 
correlation functions are of the form [2^ [3l| 

a,j,(t, t') = 2?7fe-f'(*+*') + n^7y(l - 7y)fe~l*-''l (17) 

Following the method developed in Ref. [SOj], one can 
perform a homodyne detection of the outgoing field 
fluctuations with a temporally matched local oscillator, 
^Lo{t) ^ e~^* [3^ . The measured quantity is then 
Ax^ = 1 + 77(71^ — ni^) + 7^(1 — ri)nlp, where nil/ are 
the initial and final thermal excitation numbers. When 
the cooling is optimal (77 ~ 1), one indeed measures the 
change in temperature of the mirror: Ax^ ~ 1 + — 

We have presented a general theory for the self-cooling 
of a mechanical oscillator to the ground state, which al- 
lows for deriving analytical results for the mirror final 
temperature and provides a simple interpretation of the 
cavity mediated self-cooling process. When the cavity is 
suitably detuned the initial thermal noise of the mirror 
vibrational mode is essentially transferred to the cavity 
field mode. This noise exchange reflects in the field leak- 
ing out of the cavity, providing useful information on the 
mirror temperature. 
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